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The laws of thermodynamics put limits to the efficiencies of thermal machines. Analogues of these 
laws are now established for quantum engines weakly and passively coupled to the environment 
providing a framework to find improvements to their performance. Systems whose interaction with 
the environment is activefy controffed do not faff in that framework. Here we consider systems 
activefy and focaffy coupfed to the environment, evofving with a so-caffed boundary-driven Lindbfad 
equation. Starting from a unitary description of the system pfus the environment we simuftaneousfy 
obtain the Lindbfad equation and the appropriate expressions for heat, work and entropy-production 
of the system extending the framework for the anafysis of new, and some afready proposed, quantum 
heat engines. We ilfustrate our findings in spin 1/2 chains and expfain why an XX chain coupled 
in this way to a singfe heat bath refaxes to thermodynamic-equifibrium white and XY chain does 
not. Additionafly, we show that an XX chain coupfed to a left and a right heat baths behaves as 
a quantum engine, a heater or refrigerator depending on the parameters, with efficiencies bounded 
by Carnot efficiencies. 

PACS numbers: 05.70.Ln, 05.70.-a, 03.65.Yz 75.10.Pq 


Introduction 

Considerable experimental progress in various physical systems has been achieved toward the goal of controlling the 
dynamics of open quantum systems and their interactions with the environment [H-Q . For quantum computations or 
digital coherent quantum simulations, one may wish to have a system that is well isolated from the environment. For 
dissipative variants of quantum computations Q or creating new scenarios for non-equilibrium many-body systems, 
one would need to engineer the coupling to the environment. Recently, a setting in which the quantum system of 
interest interacts at its boundaries with an external quantum probe such that their coupling can be localized and 
can be switched on and off repeatedly with a controlled and well-defined state for the probe prior to the interaction 
has been experimentally realized This repeated interaction scheme has also been theoretically studied [1, 01 . 
Importantly, the dynamics in an appropriate limit is a boundary-driven Lindblad equation. In this article, we explore 
the question of what is the thermodynamic cost of having such operations on an open quantum system and what are the 
thermodynamical quantities, such as heat and work that will determine the efficiency of quantum engines operating 
in this manner. Boundary-driven Lindblad equations have been intensively studied theoretically, particularly for 
one-dimensional quantum chains [Ml, and powerful techniques have been developed to find their non-equilibrium 
steady states (NESS) [l^ - [ll . These equations are also frequently used to describe quantum engines [20l - l^ and other 
complex open quantum systems coupled to one or several environments [23 - [T7I| because they are easy to implement. 
Nevertheless, a boundary driven Lindblad equation does not correctly describe a quantum system passively and 
weakly coupled to a heat-bath as often occurs in natural systems. It was pointed out recently that inconsistencies 
with the second law of thermodynamics may arise in this case and a careful examination of the coupling between 
a quantum refrigerator and the heat-baths [ll reveals why boundary driven models are inappropriate for these 
situations. For a system passively and weakly coupled to one or several heat-baths the master equation derived in 
the Born-Markov-secular approximation [s^ yields a proper description of the system and the correct balance of heat 
flows and irreversible entropy production. 

Thus, for our study, we consider explicitly the active (time-dependent) type of interaction between the system and 
the environment implemented in Q and the model developed in [1,0|. We apply the results of [Sll, to derive the 
appropriate thermodynamical quantities and, in particular, we focus in the limit where the system is described by a 
boundary driven Lindblad equation. 

Our main result is that driving at the boundaries, even though it looks like a work-free operation, actually might 
bring work to the system. We illustrate our findings on boundary-driven spin 1/2 chains coupled to one or two 
heat baths. We show that an XX spin 1/2 chain coupled in this way to a single heat bath relaxes to thermodynamic 
equilibrium while an XY does not because it is driven out of equilibrium by the power produced by the coupling to the 
heat bath. When two baths are connected to the chain, we observe that for different parameters, the chain operates as 
a quantum heat engine, refrigerator or heater, and we determine their efficiencies in the simple case of a chain of two 
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sites. The rest of the paper is organized as follows. We start by reviewing first the thermodynamics of Markovian open 
quantum systems in the weak coupling limit and second a formulation where the “universe”, system plus 

the environment, evolves unitarily. After that, we consider the repeated interaction scenario for the system and the 
environment from which the boundary-driven Lindblad equation and the appropriate thermodynamical quantities for 
the open system are obtained. Then we illustrate our results in XX and XY spin 1/2 chains and offer our conclusions. 
Finally we have collected in section Methods some details of the calculations. 


Thermodynamics of open quantum systems 

open system weakly and passively coupled to the environment: 

Let us briefly review the usual formulation of thermodynamics in open quantum systems . Consider an open 

system described by a master equation in the Lindblad form 

^=-^[Hs{t),ps]+Y.Vr{ps) ( 1 ) 

r 

where the environment consists of several heat-baths r whose action on the system is represented by the dissipator 

T^rip) = 

with [•,•] the commutator and {•,•} the anti-commutator. The operators L/ are system operators and represent the 
action of the environment over the system. When this equation is obtained from the weak coupling limit for a time 
independent system, one finds global Lindblad operators L/ that are eigen-operator of the Hamiltonian Hs [^ . For 
simplicity, we consider that the system can only exchange energy and no particles with the environment. 

Now consider the internal energy U{t) = ti{Hs{t)ps{t)) and the entropy S{t) = —kBTT{ps{t)\nps{t)). The first 
law of thermodynamics U = W + Q splits the rate of change of internal energy in two, power W (t) and heat flow 
Q{t) = J2rQr{t) with one contribution per heat-bath. For system passively and weakly coupled to the heat-baths, 
these quantities are defined as 


W{t) = tT[Hsit)psit)], Qrit) = iT[Hs(t)Vr{ps(t))]. (2) 

In section Methods: Heat from a given reservoir in the weak-eoupling limit we justify these definitions. Note that if the 
Hamiltonian of the system is time independent, no work can be performed on the system and only heat is exchanged 
with the baths. In that case the system will typically reach a steady state. Consider now this to be the situation. 
The second law states the positivity of the entropy production {diS/dt > 0), which is the difference between the 
time-derivative of the entropy S{t) = — ksT^'i^i'Hrips) ln/9s(t)) and the entropy flow from the environment to the 
system J^rl^rQr, 


^ = S-J2 PrQr = -kB Y, Tr(P,(ps)(lnps - \nu:p^{Hs))). (3) 

r r 

The canonical distribution ujp^(H) = exp(—appears in the last equality of Eq.([3]) due to the definition 
of heat that we plug in the first equality in Eq.([3]). The second law diS/dt > 0 in Eq.([3|) holds if for every r, 
PS = —i\Hs,Ps] +'Hr{ps) relaxes towards the unique equilibrium state ijjp^{Hs). This is the local-detailed-balance 
condition i.e. if a single heat-bath is in contact with the system detailed balance as defined in holds. This 

property of the dissipators 'Dj. is satisfied in quantum master equations obtained in the weak-coupling and with the 
Born-Markov-secular approximation (global Lindblad equation). This framework has been applied successfully to the 
study of thermodynamic properties and efficiencies of engines [2^, [39l - l^ . 

In boundary-driven systems the Lindblad operators L/ act locally on the boundaries of the system and in general the 
corresponding Lindblad equation does not satisfy local-detailed-balance. We come back to this point later. Following 
recent developments in the physics of non-equilibrium systems that have emphasized the importance of time reversal 
symmetry at the microscopic level of description , a formulation of quantum thermodynamics in which the system 
plus the environment evolves unitarily has been proposed (3lj |. We consider this framework to analyze boundary 
driven systems. 
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“universe” under unitary evolution 

Let a system and an environment with Hamiltonians Hs(t) and He (time independent), respectively, coupled by 
an interaction potential V{t) evolve with total Hamiltonian -fftot(t) = Hs{t) + He + V{t). The environment might 
consists of several heat baths pe = with = e~^^’’/Zr the initial density matrix for the reservoir 

r. Initially, the system and heat baths are uncorrelated /9tot(0) = psiO)'^P e- For arbitrary strength coupling between 
the system and environment [Slj, the internal energy is defined by E{t) = Tr(ptot(t)(F^s(t) + ^(t))); and the first law 
relates its changes to work and heat AE{t) = W{t) + Q{t) with the work W{t) = Tr(ptot(t)L^tot(t) — Ptot(0)Htot(0)) 
performed on the system in the time interval [0,t], which is also given by 

W{t) = f dt'i:,[pUt'){Hsit') + V{t')))] (4) 

Jo 

and the total heat flow Q{t) = J2r Qr{t) split in reservoir contributions 

Qr{t) = Tr{HrPtot{0)) - Tr{HrPtot{t)) ( 5 ) 

given by minus the change in energy of the r-reservoir. 

Considering S{t) = — A:BTrs(ps(t) Inps(t)) as the thermodynamic entropy of the system and AS{t) = S{t) — S'(O), 
it is found that AS{t) = AeS{t) + A^S(t) with the entropy flow AeS{t) = J2r PrQr{t) determined by the heat flows 
in Eq.([5]) and the entropy production [3l| 

AiS{t) = D{ptotit)\\ps{t) 0 Pe) > 0 , ( 6 ) 

with D{a\\b) = Tr(alna) — Tr(aln6). Unitarity, expressed through the invariance of Tr(ptot(t) Inptot(t)) under the 
time evolution of the full system, plays a crucial role in the splitting of entropy change in the entropy flow and a 
positive entropy production. In the weak-coupling limit H —>■ 0 and assuming that the open system satisfies a Lindblad 
equation obtained from the Born-Markov-secular approximation [s^, the rate of entropy production diS/dt > 0 and 
the above expressions for work and heat take the standard form given in Eg. (1^ and Eq.([2]) respectively. This is 
shown in section methods by considering the method of full-counting statistics |43| . However, the Lindblad models 
investigated in [THIOL [12 h 1M are not obtained from the weak-coupling limit and do not satisfy local-detailed- 
balance. Thus to obtain the appropriate expressions for the thermodynamical quantities in boundary driven systems 
we apply in the next section the previous formulation, in particular Eas. (|4l5l6l) . to a system plus environment evolving 
unitarily in which the reduced density matrix for the system satisfy a boundary driven Lindblad equation in an exact 
limit. 


The repeated interaction scheme 

Let us consider a finite system with time-independent Hamiltonian Hs and left (L) and right (i?) reservoirs composed 
of an infinite set of identical non-interacting finite systems with Hamiltonian H”, i.e.. Hr — where r is L 

or R. Each iL" interacts with the system for a time span t. This interaction is always of the same form, but to 
emphasize that interactions occur with different copies H^ + H"^ in different time intervals, we write it as V(t) = if 
t € [(n — 1)t, tit] with H” = -b Vr . At t = 0, the system and reservoirs are decoupled, i.e., ptot = Ps(0) 0 pe, with 
ps(0) arbitrary and pE = Pn, where Pn = a;/3j;,(7L£) 0 At t = 0+, the system begins to interact with the 

first copy -I- H^, and after a lapse of time r, the state of the total system is ptotij) = Ui [ps(0) 0 pi]Ul ® P 2 ® po---- 
Then, at t = t -|- 0, the interaction with the first copy is replaced by an interaction with the second copy for a time 
T and so on. A recursion relation for the state of the system is obtained @,0 by tracing out the nth copy of the 
environment (denoted as Tr„) 

ps{nT) = Tr„[{7„(ps((n - l)r) 0 pn)Ul]. (7) 

The unitaries are This is the repeated interaction scheme. Eor simplicity we considered 

only two heat-baths but the generalization to several reservoirs is straight forward. 

Let us consider the change of thermodynamical quantities in the time intervals of length r. Crucially, due to the 
resetting of the heat baths, the interaction term is time dependent. According to Eq.Q for time-independent Hs, 
work is performed at the discrete times nr where the interaction between the system and the environment changes 
because the copy in interaction changes. Performing the integral in Eq.Q between an initial time nr — e and a 
final time nr -b e, we obtain AWnr = Trtot([f^"~'’^ — C^Jptot) in the limit e —>■ 0. We simplify this expression with 
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the standard js^ assumption that Trr(ld"a;^(iJ^)) = 0. This condition will be repeatedly used; it allows us to split 
AWnr = AWL + AWr (we drop the index nr) with 

Aty, = -Tr(y”C/„ps([n - 1 ]t) ® pnUl). (8) 

We use Tr^. to denote the trace over the r = L or r = R system and Tr to denote the full trace. 

The heat flow from the bath to the system in the time interval of length r where the system interacts with the nth 
copy is evaluated from Eq.® 

AQ, = Tr(iJ,"[p„-p;]), (9) 

where p'^ = Tr5(17„p5([n — l]r) ® PnUX) is the density matrix of the nth copy of the environment at the end of the 
interaction with the system. 

The entropy production AiS* in the time lapse t is obtained from Eq.®, and after some manipulations [3ll,[32|, it 
can be written as the sum 

A,S = D{p'J\p^) + I{S' ■.n')>Q 

where the mutual information I{S' : n') = S{ps{nT) + S{p'^) — S'(C/„ps([n — 1 ]t) ® PnUX) quantifies the correlations 
built up between the system and the nth copy after time r. Note that D{p'^\\pn) > 0 and I {S' : n') > 0 and vanishing 
entropy production requires p'„ = pn and the absence of correlations between the system and the copy I {S' : n') = 0. 
Note that because before the interaction the state of the system is arbitrary and uncorrelated with the product of 
thermal states of the copy, the theory of [3l|, [s^ applies independently of the correlations built between the system 
and previous copies. 


Heat, work and boundary-driven Lindblad equation 

The index n is associated with the copy that interacts in the interval of time [(n — 1)t, nr], but the copies are all 
identical prior to the interaction (a tensor product of two canonical distributions) and the interaction is always 
of the same form. Because no confusion will arise, we drop the label n and denote the interaction V = 
Hamiltonian of the bath copy Hr and the state p„ = ® with = ujp^{Hr). It was shown @,0 that for Vr 

that satisfies Trj.[IAw/ 3 ^] = 0 and whose strength is scaled with r as W = '^rl\pT^ the system evolution Eq.0 in the 
limit r —>• 0 converges to a Lindblad evolution (see methods) 

PS =-i[7Ls,ps] + ^T’r(ps) (10) 

r 

with 'Dr{ps) = Trr[ur(ps 0 — ^Trr{u^,ps (8) This equation applied to particular systems provides 

boundary-driven Lindblad equations. 

Consider now Wr = AMA/r and Qr = AQrjr with AWr in Eq.® and AQr in Eq.([9]). In the limit r —0 with 
V = vjy/r, we obtain (see methods) 

Wr = Dr{Hs + Hr), Q, = -I4,(if,) (11) 

where Dr{A) = ^^[{vrAvr — \{Vr, A})ps{t) ® Note the first law '^riQr + Wr) = {Hs)t, where {Hs)t = 

Tis{Hsps{t)) = 'Y^r^r{Hs)- Finally, we express the entropy production rate as the difference between the time 
derivative of the von Neumann entropy and the entropy flow 

^ = -Ws{V{ps{t))lnps{t)) > 0 (12) 

r 

where the first term is computed using Eg. (1101) with 7) = second term is computed from Eg. dill) . 

Egs. dllll2ll provide appropriate thermodynamic expressions for systems evolving with Ea. dl0|) . Now we illustrate our 
findings in spin 1/2 chains. 


Spin models 


Consider an XY spin 1/2 chain with Hamiltonian 

N 


1 


— 2 X/ ^ {JxCTj CTj+i + 

^ i=i i=i 


7V-1 


( 13 ) 
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In the repeated interaction scheme we consider the couplings 


= Jdalal + alal), Vr = 


( 14 ) 


to a left r = L and a right r = R spin 1/2 reservoir copy with Hamiltonians Hr = hrl^lGr, and we take Hr = hi and 
hR = hjq- To obtain the boundary-driven Lindblad model, we scale Jr = iJXr/r. The canonical density matrices 
are fully characterized by the magnetization Mr = Tr(tT/a;^^) = — tanh(/3r,hr/2). 

Evaluating the second term on the right-hand side of Eq. (na yields the dissipator in the Lindblad from T’r(p) = 
7 r with 7 /^ = Ar-(1 ± Mr), Lf = af and where af = (trj ± i(T^)l2. Note 

that 

This system does not satisfy local-detailed-balance with respect to the Gibbs state, i.e. jZr is not 

the solution of 0 = —i[Hs,p] + TJr{p) with r either R or L because Vri^Pr) ^ 0- What can be shown is that these 
dissipators thermalize the single spin in the boundary if we disconnect it from the rest of the chain. Indeed let us 
consider the L dissipator 


T’l(p) ([2cr/'pcri - {cTi ut,p}] + e^^'^^[2ui pat -Wtcri ,p}]) 

upon evaluation we see that = 0. This is the generic situation in boundary driven Lindblad systems. 

The expression for power and heat Eq. m can be evaluated using the system hamiltonian Ea. (ll3L the coupling 
Ea. (fT4l) . the bath hamiltonian /iru//2 and the corresponding One obtain (we take Xr = Xr = X) 


and 


Ql = ‘2hRX{MR — TTs{alps{t))) 


(15) 


Wl = 2ATrs((J,aJ^a^ + 4afa|)ps(t)). 


(16) 


Replacing the indices {L, 1, 2} by {i?, N,N — 1} in Eas. (ll5ll6l) one has the corresponding Qr and Wr. To compute 
this quantities, we obtain psit) by solving the Lindblad equation Q. 

Consider the case in which the system interacts with one bath (for instance the left bath, but we drop the label L). 
In general, two situations can occur: the system relaxes to thermodynamic equilibrium in which all current vanishes 
or the system reaches a NESS if it is externally driven. 

XX chain coupled to one bath: An XX spin chain {Jx = Jy) in a uniform magnetic field hi = h coupled to a 
single bath relaxes to equilibrium: the entropy production rate, heat flows and power vanish. The equilibrium density 
matrix is not generally a canonical distribution but rather, as one can prove, is given by a generalized Gibbs state 
u}p{Ho) with Hq = ^ which is a conserved quantity, i.e., [Hr, Hq] = 0. This state is a product state of the 

canonical density matrices top for each spin of the chain and all equal to the one of the reservoir copy. Therefore, 
I {S' : n') = 0 and p'„ = p„, i.e., diS/dt = 0. Eigure[T] illustrates the relaxation to this equilibrium state by depicting 
the decaying power, heat flow and entropy production rate. 

XY chain coupled to a single bath: For an XY chain, we found that the system reaches a driven NESS. In this 
NESS, entropy production is strictly positive and constant, and because {Hr) = 0, the first law gives W = —Q. 
Furthermore, by combining the first and second laws, we have that [3W = diS/dt > 0 because in NESS, .5 = 0. See 
Figure [T] 

XX chain coupled to two baths: Consider a hot left and a cold right heat baths (/ 3 l < (^r) connected by an XX spin 
1/2 chain with the Hamiltonian in Ea. (fT51) with Jx = Jy = J■ The NESS in the special case of a uniform magnetic 
field was analyzed in Q. The power and heat from the reservoir to the system are given by Eas. (|15ll6D . In Figure [21 
we plot Qr, Qr, W = Wr + Wr and diS/dt in the NESS as functions of hR. We observe that the heat flows can 
change sign and that for hR = hR, they have opposite signs, i.e., Qr + Qr = 0, which means that W = 0. We also 
observe in Figure H] that diS/dt > 0 and vanishes only when ^RhR = (dRhR, that is, the second law holds even when 
heat flows from cold to hot, as is the case for hR > fdRhR/{H r, a situation that would appear to be a contradiction to 
the Clausius statement of the second law if we do not realize the presence of W. 

The previous numerical study of boundary-driven spin chains can be complemented with exact results for power 
and heat in a two-site boundary-driven spin chain obtained from a full analytical solution of the NESS (see methods). 
In the NESS, the expression for power Eq. (HU) and heat Eg. (1151) can be written in terms of the spin current (3| 

, _ U\Mr-Mr) 

{hR - /lfl)2 + 16J2 + 16A2 
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FIG. 1: As a function of time t plots of W (blue), —Q (black) and diS/dt (red) for an XX (Jx = 1 = Jy) and —W (blue, 
dashed), +Q (back, dashed) and —diS/dt (red, dashed) for an XY {Jx = 1 — 0.5Jy) chain. In both cases, the chain has N = 5 
sites with /li = 1, i = 1, 5 coupled with A = 1 to a single left bath of /I = 1 and h — 1. 


as Ql = —hLjs, Qr = hpijs and W = {h^ — hfi)js. Thus, for = h/j, there is no power, but as the previous 
expression shows, this does not mean that the spin current vanishes. Moreover, the entropy production rate in the 
NESS is 

7 Q 

= {PLhL - PRh[i)js 

i.e., the spin current js and the affinity {Pl^l — ^RhR) characterize the rate of entropy production in the NESS, and 
because Mr = — tanh(/3r/ir/2), the sign of the entropy production rate is given by {^LhL — /3_R/i_R)(tanh(/3L/iL/2) — 
tanh(/3_R/ii{/2)) > 0, where the equality holds only if = PrUr. Let us end this analysis by noting that for 

/3l < this system behaves as a heat engine for Pl/Pr < hR/hR < 1 with efficiency rj = —WIQl = 1 — hR/hR < 
1 - Pl/Pr = VC, as a refrigerator for hR/hR < Pr/Pr < 1 with efficiency = QrIW = IjihRlhR - 1) < 
^UPlIPr ~ 1) = ^ heater for Hr/Hr > 1. Note that the efficiencies are independent of temperature. These 

are steady-state operating engines analogous to those in (4H| . 



FIG. 2: For a N = 5 site XX chain with Jx = Jy = 3, h 2 = /13 = /14 = 5, = hn = 2, Pr = 0.8, Pr — 1.2, and A = 1, we depict 

li/^®®® and (diS/dt)^®®® as a fnnction of Hr = hi. There are two special values for Hr. At Hr = 3, where 
PRhR = PrUr, all quantities vanish (eqnilibrinm state). At Hr = hn = 2, Q^®®® = —and thus IF”®®® = 0 (non-driven 
steady state). 


Discussion 

In conclusion, the repeated interaction scheme provides a physical description of a system interacting with an 
environment that, in an appropriate limit, provides a boundary-driven Lindblad equation for the system. The Lindblad 
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operators that appear in this equation are determined by the interaction of the system with the environment, the 
Hamiltonian of the copies that form the bath and, importantly, by the fact that it is refreshed constantly. By computing 
the thermodynamical quantities for the full system plus the environment, one can derive the corresponding expressions 
for the boundary-driven model. One important observation is that due to the refreshing of the reservoir, work is done 
or extracted by the external agent in charge of this refreshing. This power drives the system out of equilibrium. Note 
that this power appears even if the system Hamiltonian and Lindblad operators are time independent. We applied 
our results to spin chains. In the single bath case, we found that an XX spin chain with a homogeneous magnetic 
field relaxes to thermal equilibrium, i.e., a state with zero entropy production, while an XY spin chain reaches a 
driven NESS, a state with a non-zero entropy production diS/dt = fiW > 0. In the two heat bath case, the XX chain 
for different temperatures fin ^ (3^ and a homogeneous magnetic field reaches a non-driven W = 0 NESS and an 
equilibrium state for = pRhn where the entropy production rate, power, heat flows and spin currents vanish. 

For inhomogeneous magnetic fields, the chain reaches a driven W 0 NESS. Jumping to a broader context, this 
work shows that the knowledge of a Lindblad equation for an open system does not determine the heat flows or 
other thermodynamical quantities. These quantities also depend on the properties of the environment and how the 
system is coupled to it. Here, we have obtained appropriate expressions for heat flows and power for interactions with 
an environment of a type recently implemented in a laboratory Q. But when the reservoir is weakly and passively 
coupled to the system, i.e. there is no work cost in achieving the coupling, the system is appropriately described by a 
global Lindblad equation and the thermodynamical quantities by Eq. ©■ Finally, this work is also an extension 
of quantum thermodynamics to a class of open quantum systems without local-detailed-balance. 


Methods 

We provide here some details of the calculations mentioned in the main text. 


Work, heat and boundary-driven Lindblad equation from the repeated interaction scheme 

For completeness we derive Eq. (US and Eq. (HU of the main text. Consider Apsinr) = ps{nT J- t) — psinr). We 
have from Eq. ([7]) of the main text that 

Aps(nr) = TrE[Ups{nT) 0 peU'^ - psinr) (g) pe] (17) 

where we have dropped the label n from {7„ and in Eq. © because the copies are identical and the interaction = 
K" i® always of the same form. The trace Tr„ over the state Pn = Pe = ® ^Pr i® denoted Trg. The unitary 

U = jjj (|T71) is expanded for small r considering the scaling V = vjy/r and Hq = Hs + Hr -|- Hr 

U = I - ivT^ - (iHo + y)t - i {Ho, n} r® -k O(r^). (18) 

Now, because Tt:e[vpe\ = 0 the leading order in the right hand side of (IT71) is 0{t). Thus, we divide by r and take 
the limit r —>■ 0 and n —> oo such that t = nr and obtain 

PS = -ilHs: Ps] + Tri5[n(ps (g) pe)v] - ^Trsl'u^, ps ® Pe} 

where the equality Tve{[He, Pe\) = 0 was used. 

Now we use Tr^ to denote the trace over the r = L oi r = R system and Tr the full trace. Because v = Vr and 
Tir[vrUJp^] = 0, it is possible to split the last two terms in contributions for each reservoir giving Ea. dTOl) in the text: 

PS = -i[Hs,ps]+'^Vrips) (19) 

r 

with Vrips) = Tirlvrips <8 ujp^)vr] - lTrr{n^,ps (g W/J,.}- 

We continue with the derivation of Eg. lfTTI) of the main text. Let us start from AQr = Tr(iJ"[p„ — p(i]), i.e. Eq.(j9l), 
where p'^ = Trs(C7„ps([n — 1 ]t) ig) pnU^). Dropping as before the label n, in the limit V —>■ vj^Jr and r —>■ 0 we can 
replace U by (HHl). The leading order of AQr is 0{t) 


r - ^{Vr.Hr})psi[n-l]T) < 


' Pe 


AQr = —rTr I (VrHrV, 


(20) 


Qr = -Tr (^VrHrVr - ]^{vl, Hr})ps{t) 0 

Consider Eq.dH) now i.e. AWr = —Tv{VrUnPs{[n— 1]t) (g) PnUX)- As before we drop the label n. The leading order 
is also 0(t) but we need U up to 0(r^/^)because V is 0(l/v^), pe) = 0 and ^ie{[{Hs + He),V]ps^ ps) = 0. 

We obtain 


AWr = rTr (vriHs + Hr)t 


- 2 ^Vr, Hs + Hr})ps{[n - l]r) 0 W/j, 


or Wr = AWrjr 


Wr = Tr 


^^Vr{Hg Hr)Vr 




Hs + Hr})ps{t)(S)UJi3^ 


Expressions (uni) and m correspond to those in Eq. m from the main text. 


( 21 ) 


The two spin XX chain with inhomogeneous magnetic field 

Consider a XX two sites spin chain and the corresponding Lindblad dynamics Eq. © with Hs given by Eq. m 
main text (with = Jy = J, hi = hs and /12 = hn) and the Lindblad dissipator 

t^Ap)= Y. li^mpLA -{lAl>},p}] 

mg{+,—} 

with jA = '^'■(l ± Mr), = af and = af where = (crj ± ia^)l2. This system is fully characterized by the 
correlation functions {(T})t, (uDt, Y = j((Tfcrf — (j}(J 2 )t and X = (afa^ + where (■}t = Trs(-ps(t))- They 

satisfy a close system of equations: 

^ = -4AX - i(hn - hL)V (22) 

at 


= 4A(Mi - {a})t) + 2iJY 


(23) 


= 4X{Mn - (aDt) - 2iJY 


(24) 


— = i{hL - hR)X - - K)t) - 4Ar (25) 

From Eas. lfTKl and fTHll in the main text we note that Wr = 2XJX = Wr, while the first term in the right hand side 
of (|23|) is 2QL/hL and the corresponding term in (l24l) is 2QR/hR. Moreover the spin current Q is js = iJY. In the 
steady state the left-hand-side of the system (|22I23I24I25I) vanishes and W = Wr -\- Wr = [hs — hR)js, Ql = —hsjs 
and Qr = hRjg. The current given in the main text is obtained by solving the full system in the NESS. 


Heat from a given reservoir in the weak-coupling limit 

Consider a system coupled to several reservoirs as discussed in “universe” under unitary evolution. The heat that 
comes from one of them, for instance the r = L reservoir is Ql = Tr[iLL(ptot(0) — ptot(t))]. The methods developed 
in full counting statistics gives Ql = *9aG(A)|a=o = —iYT{td\CxPs{t))\>^=o where £\p = —i{HAtP — pHiJ}) is a 
modified evolution super-operator with iLtot —t H^t = ^ When this modification is done for 
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a system in the weak coupling Born-Markov-secular approximation one obtain [4^ C\p = —i[Hs,p\ + 
where only the dissipator associated to the r = L reservoir depends on A as 

Vl.xY = EE - {A^,Y}) + {2A^YA\‘^ - {A\‘^A'^,Y}). (26) 

I UJ 

Here are system eigen-operators obtained from the coupling of the system to the left reservoir [s^, |4^ and 
h'l~{uj) = {uj). A slow time dependence of the system can be included, see [d^. From Ea. d^ we obtain 

- z(5ApL,A)U=oi" = {hi{uj)A\^YAt - hr{oj)ArYAl‘^) . (27) 

Thus Ql = —*Tr[(c}AT’L,A)ps(i)]|A=o where we used that in this limit the dynamics is Markovian. We have to 
compare this with the heat flow defined in section “open system weakly and passively coupled to the environment”, 
Qrit) = tv{Hs'DL{ps)) = It:{ps'D\{Hs)), where the dissipator Vl in the same weak coupling Born-Markov-secular 
approximation is given by I7 l,a= 0) from which we compute 

v{{Hs) = 2^cc {hi{uj)A^A\“^ - hY{oj)A\“^At) . (28) 

l,UJ 

To obtain this we used [Hs,A‘^^A‘p] = 0 js^. Taking the trace in Eg. (071) and in Eg. (051) the desired equality 
—iTr[(9AT>i_A)ps(i)]|A=o = ^^{HsT^LiPs)) is found. Now, since the heat flow to a system weakly and passively 
coupled to the L heat-bath is given by Qrit) = tr(i7sT>L(ps)), the corresponding definition for work follows and the 
entropy production given in Eq.(I3]) as well. 
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